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We study the impact and subsequent retraction dynamics of liquid droplets upon high- 
speed impact on hydrophobic surfaces. Performing extensive experiments, we show that 
the drop retraction rate is a material constant and does not depend on the impact 
velocity. We show that when increasing the Ohnesorge number, Oh = rj/y/pRij, the 
retraction, i.e. dewetting, dynamics crosses over from a capillaro-inertial regime to a 
capillaro-viscous regime. We rationalize the experimental observations by a simple but 
robust semi-quantitative model for the solid-liquid contact line dynamics inspired by the 
standard theories for thin film dewetting. 



1. Introduction: Drop Impact on Solid Surfaces 

Drops impacting onto solid surfaces are important for a large number of applications: 
for instance, almost all spray coating and deposition processes rely ultimately on the 
interaction of a droplet with a surface. A large variety of phenomena can be present 
during drop imp acts, from sp lashes to spreading, and from large wave surface deformation 
to rebound (see l)Reinlll99^) and references therein). 

Research on drop impacts has a long history, starting with the p ioneering studies of 
Worthington and later on with the famous photographs of Edgerton l)Worthingtorjll~876t 
lEdgerton fc Killianlll954 - Most of the previous work on drop impact focus ed on deter- 
mining the maximum diameter a drop is cap able of covering upon impact l|Fukai et all 
ll993tlRoisman et alJl2002HClanet et a/.ll2004|) . However, the practical problem of deposi- 
tion can be very different if one wants to efficiently deposit some material on the surface. 
This is especially grave when the surface is not wetted by the liquid, as is illustrated 
by the high-speed video pictures in Fig^for the impact of a water droplet. It can be 
observed that the drop expands rapidly, due to the large speed with which it arrives at 
the surface. However, due to the hydrophobicity of the surface, subsequently the drop 
retracts violently, leading to the ejection of part of the droplet from the surface: we ob- 
serve droplet rebound. It is this "rebound" that is the limiting factor for deposition in 
many applic ations, for instance fo r the deposition of pesticide solutions on hydrophobic 
plant leaves ^Bergeron et aZ.ll2000^l . We study here the impact and subsequent retraction 
of aqueous drops onto a hydrophobic surface, and seek to understand the dynamics of 
expansion and retraction of the droplets. 
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Figure 1. Temporal evolution of the contact radius of droplets upon impact and retraction. 
The radii are normalized by those of the spherical droplets before impact. The pictures show 
the shape of the droplets at the different stages of retraction. Droplet radius is 1 mm, impact 
speed is 2 m • s _1 : a) pure water, b) viscous water-glycerol mixture, viscosity 50 mPa • s. 

In general, these problems are difficult because for most practical and laboratory sit- 
uations, three forces play an important role: the capillarity and viscous forces, and the 
inertia of the droplets. We try and disentangle the effects of the three forces here by 
performing systematic experiments, varying both the importance of viscous and inertial 
forces. 

We provide experimental evidence for the existence of two distinct retraction regimes. 
In both regimes, capillary forces are the motor behind the droplet retraction, and are, 
for the first regime countered by inertial forces. In the second regime the main force 
slowing down the retraction is viscous. We also show that, perhaps surprisingly, the drop 
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retraction rate (the retraction speed divided by the maximum radius) does not depend on 
the impact velocity for strong enough impacts. The dimensionless number that governs 
the retraction rate is found to be the Ohnesorge number, Oh = r\j y/pRij, with 77 the 
viscosity, p the liquid density, R\ the impacting drop radius, and 7 the surface tension. 
The Ohnesorge number therefore compares the dissipative (viscous) forces to the non- 
dissipative (capillary and inertial) forces. The crossover between the two regimes is found 
to happen at a critical Ohnesorge number on the order of 0.05 . 

In order to develop a better understanding for the different regimes that are encoun- 
tered, particularly the retraction dynamics in these regimes, we propose two simple hy- 
drodynamic models inspired by the standard description of thin film dcwctting dynamics. 
These simple models provide a simple but quite robust picture that allows us to ratio- 
nalize the retraction rate in both regimes. 

In order to be able to say something about the speed of retraction, one also needs to 
understand the maxim um radius to which the droplet expands. Combining our results 
with those obtained by l|Clanet et a/.ll20o"^ for the maximum radius, we propose a phase 
diagram delimiting four regions for the spreading and retraction dynamics of impacting 
drops. 

2. Drop retraction dynamics: Generic Features 

As the impact dynamics of liquid droplets on a solid surface happens usually in a few 
tens of milliseconds, we use a high-speed video system (1000 frames/second, Photonetics) 
to analyze the drop-impact events. When necessary, we use an ultrahigh-speed system 
allowing to go up to 120,000 frames/second (Phantom V7). We study aqueous drops 
impacting on a solid surface; the surface we used is Parafilm, which provides us with a 
hydrophobic surface (receding contact angle for water #r ~ 80°). In addition, the surface 
has a low contact angle hysteresis with water, and allows us to obtain highly reproducible 
results. The liquids we used are different water-glycerol mixtures. Varying the glycerol 
concentration, we vary the liquid viscosity, keeping the liquid density and its surface 
tension almost constant. For the highest concentration of glycerol, the surface tension 
has decreased from 72 (pure water) to 59 mNm^ 1 , whereas the density has increased 
to 1150kg/m 3 . The viscosity is varied between 1 and 205 mPas. Viscosity, density and 
surface tension were measured before each impact experiment. Drops were produced 
using precision needles, and the initial radius of the drops Ri have been systematically 
measured on the images (1.1 < R\ < 1.4 mm). From the high-speed images such as the 
ones shown in FigQ] we follow the contact radius R in time. This section summarizes 
the results of more than 80 different drop impact experiments, each of which have been 
repeated at least two times. 

Two series of experiments were performed: first, letting the droplets fall from a fixed 
height, but increasing the viscosity, we increase the Ohnesorge number while keeping the 
inertial forces constant. The second series of experiments is performed at fixed viscosity 
and upon increasing the height from which the droplets falls; the droplet turns out to be 
in free fall (as is verified in the experiment to within a few percent) and so the relation 
between fall height h and impact velocity is simply V\ — \fgh, with g the gravitational 
acceleration. Increasing the impact velocity increases the Weber number, keeping the 
Ohnesorge number fixed, where the Weber number, We, compares the inertial forces to 
the capillary forces, We = pR\ Vf/7. 

In all that follows, we restrain ourselves to high-speed impact conditions. More pre- 
cisely, the Weber and Reynolds numbers are chosen so that We > 10 and Re > 10, where 
Re = pR\ Vi/n is the Reynolds number. This implies that inertial forces are at least one 




Figure 2. Temporal evolution of the contact radius for a water-glycerol drop Oh — 9.1 10~ , 
Ri = 1.2 mm (a) contact radius vs. time, (b) contact radius normalized by the maximum spread- 
ing radius vs. time Impact velocities : X: Vi = 2.4ms -1 , +: V\ = 2.2ms -1 , o: Vj = 1.9ms -1 , □: 
Vi = 1.7ms -1 , A: Vi = 1.4ms -1 , o: Vi = 1ms -1 



order of magnitude larger than both the capillary and the viscous forces. Such conditions 
imply large deformations of the drop when the liquid impinges on the solid substrate. 
On the other hand, we also restrain our experiments to impact speeds that are far from 
the 'splashing' regime in which the drop disi ntegrates after impact to form a collection 
of much smaller droplets l|Mundo et al,l\l99^tf . 

The pictures in Fig^show that two distinctly different regimes exist for the shape of 
the droplets after impact. For low fluid viscosity, we typically obtain the images shown 
in Fig^a). At the onset of retraction, almost all of the fluid is contained in a donut- 
shaped rim, with only a thin film of liquid in the center. On the other hand, for high 
viscosities the deformation of the drop is less important, and the pancake-shaped droplet 
of Fig^b) results. These visual observations allow to distinguish the capillary-inertial 
and the capillary-viscous regimes that are described in detail below directly. 

2.1. Drop Retraction Rate: influence of fall height and viscosity 

Fig. [21 summarizes the most important findings of this study. The temporal evolution of 
the drop contact radius R(t) for different impact velocities, shown in (a), is normalized in 
(b) by its maximal value at the end of the spreading i? max . Two important observations 
are made, (i) A well defined retraction velocity V lct can be extracted from each experi- 
ment; this is a non-trivial observation that will be rationalized below, (ii) Independently 
of the impact speed, all the R(t)/ R mayi curves collapse onto a single curve for different 
impact velocities. This shows that the retraction rate, rather than the retraction speed is 
the natural quantity to consider, and that this rate is independent of the impact velocity. 
These results hold for all the viscosities tested in our experiments. 
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Figure 3. Retraction Rate plotted versus Impact Weber number for various water glycerol 
droplets. X: Oh = 2.510" 3 , +: Oh = 3.910" 3 , o: Oh = 1.510 -2 , A: Oh = 1.610" 2 , □: 
Oh = 2.310" 2 , o: Oh = 7.110" 2 



In Fig. 01 we have plotted the retraction rate e = ^ rGt /i? m ax versus the impact Weber 
number, where V~ rot is defined by V rc t = max [— R(t)]. Clearly, the drop retraction rate 
does not depend on the impact velocity. One might think that the explanation for this 
observation is rather obvious: the initial kinetic energy of the droplet is transformed into 
surface energy (which fixes i? m ax/-R/ oc We 1 / 2 ), and is then transformed back into kinetic 
energy (which in turn fixes V Te t oc Vj ). This naive explanation is unfortunately wrong fro 
the following reasons. First, it has been observed recently tha t, at the onset of r etrac- 
tion, low viscosity liquids undergo vortical motion in the drop l|Clanet et ~d\ \2004 . This 
residual flow in the drop reveals that a part of the initial kinetic energy is still available 
then, and thus that a simple energy balance argument cannot work. This was indeed 
already suggested by previous o bservations of a clear disagreement between experiments 
and t he R ma , x /Ri oc We 1 / 2 law l|Fukai et all\l99?i iRoisman et oiJl2002l lOkumura et al\ 
120031) . The second reason why the simple energy-balance argument does not work follows 
directly from Fig. |3 where it is shown that the retraction rate depends on the viscosity 
and consequently that the previous inviscid picture is not correct. 

We therefore performed experiments that elucidate the role of the viscosity, or, equiva- 
lently, of the Ohncsorgc number. For what follows, it is convenient to define two intrinsic 
time scales for the droplet: a viscous one and an inertial one. The viscous time is the 
relaxation time of a large-scale deformation of a viscous drop: r v = (rjRi) /y, whereas the 
inertial time scale: t\ = {^'KpR\/ r )) 1 l' 2 corresponds to the capillary oscillation period of 
a perturbed inviscid droplet. Since rj is independent of V\ and 77, this quantity is almost 
constant for all tested drops. 

Fig. 01 shows the retraction rate, made dimensionless using the inertial time, as a func- 
tion of the Ohnesorge number. It can be observed in the figure that two different regimes 
exist for the retraction rate. The first region where the retraction rate e is independent of 
the viscosity points to an inertial regime and e oc t^ 1 . The retraction rate is consequently 
found not to depend on the impact speed, a result similar to that obtained recently by 
ijRichard et allEw^l who show that the contact time is independent of the impact speed. 
For higher viscosities, typically Oh > 0.05, the retraction rate decreases strongly. In this 
regime, capillary and viscous forces govern the dynamics: we find e oc t v _1 . 
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Figure 4. Circles: Normalized retraction rate erj plotted versus the Ohnesorge number, exper- 
imental values. Error bars represent the maximum deviation from the mean value. Full line: 
(left) eri evaluated using Eq. 13.21 (right) erj evaluated using Eq. 13.61 Dashed line :(left) Fit 
obtained taking the mean value of the five first experimental points, (right) Best fit according to 
the predicted 1/ Oh power law. 

3. Two simple models for the drop retraction dynamics 

We have consequently established the existence of two different regimes for the retrac- 
tion rate: a viscous one and an inertial one. We now develop some simple arguments 
allowing for a semi-quantitative description of the dynamics, using ideas already existing 
for the dynamics of dewetting, a problem closely related to the current one. 

3.1. Inertial regime 

We e mploy a Taylor-Culick approach commonly u sed for the inertial dewetting of thin 
films l)Tavlodll959HCulici3ll96filBuguin et a^Jll99^) to describe the drop retraction rate. 
For high-velocity drop impacts, liquid spreads out into a thin film of thickness h and 
radius i? max - The liquid subsequently dewets rapidly the surface, and in doing so forms 
a rim that collects the liquid that is initially stored in the film. The shape of the drop 
surface shape is therefore never in a steady state and consists of a liquid film formed 
during the spreading stage and a receding rim. The contact angle at the outer side of the 
rim is taken to be very close to the rec eding contact angle since viscous effects can 
be neglected small llBuguin et aZ.lll99^) . The dynamics is consequently determined by a 
competition between capillary tension coming from the thin film and the inertia of the 
rim. If we write down momentum conservation for the liquid rim: 



with m the mass of the liquid rim and Fq the capillary force acting on it, Fc ~ 
2-K~fR{t) [1 — cos(0r)]. The stationary solution of Eq 13. II can be obtained writing m(t) = 




(3.1) 
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2irpRVieth, and gives: V re t = y/lW — cos (#r)]/ (ph). Using volume conservation, h 
^RjR^naxi it follows that: 

V rrt 



Rr, 



T; 



y/n [1-cos0r)] (3.2) 



Which is the final result. Comparing with the experimental data, it turns out that this 
equation not only gives the correct scaling behavior for the retraction in this regime 
rate but also provides a rather accurate estimate of the numerical prefactor (see FigUJ. 
Indeed, the ratio between the experimental and the predicted numerical prefactors is 
found to be 0.6 Repeating the experiment for water on a polycarbonate surface, which 
changes the contact angle value to 60°, we retrieve exactly the same ratio of 0.6. 

3.2. Viscous regime 

In the opposite limit of very viscous liquids, the drops adopt pancake shapes upon im- 
pact. During the first stages of retraction, the pancake shape rapidly relaxes towards a 
roughly spherical cap, and the drop shape remains like this during the retraction since 
the capillary number is small. During the retraction, it is only the contact angle that 
varies slowly: it is mainly this slow contact angle dynamics that dictates the drop evolu- 
tion during the retraction. Contrary to the previous analysis, the slow receding velocity 
allows to assume a quasi-static dynamics for the surface shape during the retraction. 
In this regime, it is then natural to assume that the work done by the capillary force 
Fq is dissipated through viscous flow near the contact line. Since we focus our study 
on high-speed impacts, -R max is always much larger that R\ which justifies a small 8(t) 
approximation at the onset of retraction. The visco us effects near th e contact line then 
lead to the well-known linear force- velocity relation ljDeGenneslll985|) : 

F v = -^hi(j^R{t)R(t) (3.3) 

where A and A are respectively a macroscopic and a microscopic cutoff lengths. A is 
typically of the same order as the drop size ~ 1mm. A is a m icroscopic length, and is 
usually taken to be on the order of A ~ lnm ijDeGenneslllflS^ . On the other hand, the 
capillary force drives the retraction. Near the contact line it can be written: 

F c = 2tt R{t)~/ [cos 6{t) - cos R ] (3.4) 

Volume conservation gives: |7ri?j ~ j8(t)R 3 (t), where we have taken the small angle 
limit. Eg s . 13 . 31 and 13 . 41 together with the volume constraint leads to the following relation 
for the variation of the contact radius: 

R{tl [1- l0 2 (f) -CQ S (fl fl )] 9(t)^ 

R(t) (144) 1 /3l n (A/A) Tv 1 ' 

the above equation is obtained in the small angle limit and is only valid for short time 
after the onset of retraction. We estimate the retraction rate e as the maximum value of 
R(t)/R(t) so that: 

^rct /3\ 1/3 (l-co S fl R ) 5 / 3 _ x 

i?max ~ \25j 51n(A/A) Tv 1 ' ' 

Comparing again to the experiments, good agreement is found: the retraction rate is 
solely set by the viscous relaxation time r v and consequently er; cx Oh -1 . Beyond this 
correct scaling prediction, Eq. 13.61 provides a quite accurate estimate for the numerical 
prefactor as is shown in Fig 0] Indeed, the ratio between the experimental and the 




Figure 5. (a) Normalized maximum spreading radius plotted vs. the impact number, (b) -R max 
(normalized by the radius before impact) plotted vs. Weber number for small values of the im- 
pact number. Full line: power-law fit. (c) 7? ma x (normalized by the radius before impact) plotted 
vs. Reynolds number for large values of the impact number. Full line: predicted power-law de- 
pendence with power 0.2. *: 77 = W^Pa.s, +: 9.51CT 2 , o: 77 = 4.81(T 2 Pa.s, A: 77 = 2.810~ 2 Pa.s, 
□ : 77 = 10~ 2 Pa.s 



predicted numerical prefactors is found to be 1.5. Again, repeating the experiment on a 
polycarbonate surface, this ratio changes only slightly from 1.5 to 1.8. 



4. Conclusions and perspectives 

Our experiments reveal that the retraction rate is independent of the impact speed. 
To account for the retraction speed, the maximum radius to which the droplet expands, 
has to be known also. A number of studies have been devoted to the understanding of 
the maximum spre ading radius (see for instance l)Fukai et alJll993t iRoisman et aZ.ll20o3 
IClanet et al\\200^i ). However, no clear and unified picture emerges from previous ex- 
perimental investigations. A recent experimental study of i? max , combine d with recent 
theore tical ideas in the same spirit of the ones presented here was done bv lClanet et all 
l)2004|) . They obtain a zeroth order (asymptotic) description of the spreading stage, com- 
pare it with experiments and suggest that two asymptotic regimes exist for i? max . The 
first is given by a subtle competition between the inertia of the droplet and the capillary 
forces; if only these two are important, it follows that i? max /i?i oc We 1 / 4 . In the second 
regime, R max is given by a balance between inertia and viscous dissipation in the expand- 
ing droplet, leading to i? max /-Ri oc Re 1 / 5 . Consequently, a single dimensionless number 
is defined that discriminates between the two regimes: P = WeRe~ 4 / 5 referred to as the 
Impact number. The crossover between the two regimes happens at a P of order unity. 

Our experimental data are in qualitative agreement with their prediction, as is shown in 
Fig.|SJa. At low P, the scaling i? max /i?i ~ We 1 / 4 is clearly observed. However, for impacts 
corresponding to P > 1, we observe only a very slow variation of the maximum spreading 
radius as a function of P. Therefore, the relation between i? max and the Reynolds number 
is not very clear from our data (Fig [S] c) . Although the main trend is not in strong 
contradiction with the prediction i? max /i?i oc Re 1 / 5 , a power-law fit of our data gives 
exponents that are always smaller than the predicted value of 0.2. Perhaps even more 
important- in view of the small range of the maximal expansion Rxnax. t licit W6 cover- is 
that the different water-glycerol mixtures do not appear to collapse on a single master 
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curve, as would be predicted by the above argument. However, since the maximum value 
of P that we reach is on the order of 10, it may be that we have not reached the 
purely viscous regime. In that case, the capillary, inertial and viscous forces are still 
of comparable amplitude and have to be taken into account to gether. Note also that 
the more sophisticated models reviewed in lUkiwe fc Kw ok ( 200^) do not provide better 
agreement with our experimental measurements. 

Despite this small problem, we are now able to develop a simple unified picture for drop 
impact dynamics accounting for both the spreading and the retraction dynamics. The two 
natural dimensionless numbers that have been identified arc the impact number P, that 
quantifies the spreading out of the droplet, and the Ohnesorge number Oh that quantifies 
the retraction. We can thus construct a phase diagram in the experimentally explored 
(Oh, We) plane, which is shown on Fig.0 The experimentally accessible plane is divided 
in four parts, where the main mechanisms at work during the impact process are different. 
These four parts are separated by the curves Oh = 0.05 and We = Oh~ 4 ^ 3 . They are 
labeled as follows: ICCI the drop dynamics is given by a competition between inertia and 
capillarity both for the spreading and the retraction. IVCV: inertia and viscous forces 
dominate the spreading, capillary and viscous forces dominate the retraction. These two 
regimes have been studied in detail here. The two more intriguing regions are the IVCI 
(viscous spreading, inertial retraction) and ICCV (capillary spreading, viscous retraction) 
that are unfortunately difficult to explore in detail. For the IVCI- regime, the large inertia 
at impact, combined with a small surface tension, will make the droplets undergo large 
non-axisymetric deformations and they will eventually splash and disintegrate. On the 
other end of the phase diagram, the ICCV region corresponds to very low impact speeds 
and important capillary forces, implying very small deformations of the droplets. If the 
deformations are small, pinning of the contact line of the droplets will become important, 
and all our simple scaling arguments for both the maximum radius and the retraction 
rate are invalidated. 

A numerical investigation of droplet impact would be very helpful for two reasons. 
First, numerics would allow to vary R\ while keeping all the other physical parameters 
constant. This would allow to check the robustness of our results, since experimentally 
it is not easy to vary R\ over a wide range. Second, as emphasized above, the viscous 
regime for the maximum radius is difficult to characterize precisely due to the smallness 
of the variation of R max for viscous drops. If precise numerical simulations could be done, 
these different remaining problems could be resolved. 

In sum, we have studied the retraction dynamics of liquid droplets upon high-speed 
impact on non-wetting solid surfaces. Perhaps the strongest conclusion from our inves- 
tigation is that the rate of retraction of the droplet is a drop constant which does not 
depend on the impact velocity. Two regimes for the retraction rate have been iden- 
tified: a viscous regime and an inertial regime. We have in addition shown here that 
simple hydrodynamic arguments can be formulated that give very reasonable agreement 
with experiments in the two different regimes. Acknowledgments: Benjamin Hclnann- 
Moussa is acknowledged for help with the experiments. Denis Bartolo is indebted to the 
CNRS for providing a post-doctoral fellowship. LPS de TENS is UMR 8550 of the CNRS, 
associated with the universities Paris 6 and Paris 7. 
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